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Abstract 

We study the nonlinear transport properties of NS (normal-superconductor) 
and NSN structures by means of a self-consistent microscopic description. A 
nonzero superfluid velocity causes the various quasiparticle channels within 
S to open at different voltages. The gap reduction is very sensitive to the 
details of quasiparticle scattering. At low temperatures, superconductivity, 
some times in a peculiar gapless form, may survive up to voltages much higher 
that ksTc/e. The minimum voltage for quasiparticle transmission is shown to 
decrease strongly with temperature and with the transmittivity of the barrier. 

PACS numbers: 74.40. +k, 74.50.+r, 74.90. +n 
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During the last few years, there has been a renewed interest in the transport properties 
of structures that involve both normal and superconducting elements. This research has 
been largely motivated by the study of mesoscopic superconductivity, whose characteristic 
feature is the phase-coherent propagation of quasiparticles. Due to its fundamental charac- 
ter, the normal-superconductor (NS) interface has been the object of especial attention, and 
a number of transport anomalies have been associated to the existence of phase-coherent 
Andreev reflection.! These investigations have been accompanied by a revived interest in 
related questions such as the role of self-consistency in microscopic descriptions based on 
the resolution of the Bogoliubov - de Gennes equations,! 
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Here, Hq is the one-electron Hamiltonian, A is the gap function, and (u n , v n ) and e n are, 
respectively, the normalized wave function components and the energy of the quasiparticle 
n. Self-consistency requires 
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where g is the coupling constant and {f n } are the occupation probabilities. The implemen- 
tation of self-consistency is especially important in transport studies, since it guarantees the 
conservation of electric current. In a perfect one-dimensional superconductor, uniform 
current- carrying solutions are of the formjl 



A(x) = \A\e i2qx . 



(3) 



In equilibrium, |A| is a decreasing function of q, with a critical current above which 
|A| =oB&i Self-consistency is also essential to understand the crossover from the Josephson 
effect to bulk superconducting flow.! A smooth crossover is possible thanks to the exis- 
tence of uniform and solitonic (with a local gap depression) configurations SB For scattering 
structures, the use of scattering quasiparticle channels compatible with (3) is sufficient to 
guarantee asymptotic self-consistency.ll! 



The studies mentioned above either neglect self-consistency or assume equilibrium of 
quasiparticles. These may be reasonable assumptions in specific contexts, namely, at low 
voltages (assumption of equilibrium) or for currents much smaller than the critical current of 
the perfect superconductor (neglect of self-consistency). However, a theoretical description 
that includes both self-consistency and nonequilibrium effects is mandatory in transport 
problems involving high applied voltages and large superconducting current densities. In 
this Letter, we make a first step in this direction. First, we generalize the classic scattering 
study by Blonder, Tinkham, and KlapwijkS on the NS interface by including the effect 
of self-consistency. We will see that this inclusion introduces major changes in the I-V 
characteristics, but we will also argue that the most spectacular of them are not observable 
in practice because of the ultimate presence of a second normal lead in any real experiment .0 
For this reason, we extend our analysis of the NS interface to the study of a NSN structure. 
We compound the two interfaces by assuming incoherent multiple scattering by the two 
barriers. This simplification has the advantage that it permits us to focus on the properties 
arising specifically from the combination of self-consistency and the superposition of two 
interfaces. In this sense, we expect our predictions to be even more robust than those we 
would have inferred from coherent scattering descriptions, since these might be dependent 
on fine, not easily reproducible details of the scattering process. Therefore, the realization 
of the physics discussed in this paper does not require extremely small structures nor very 
low temperatures. The only essential requirement is that, due to the application of high 
voltages, and in the absence of dilution effects, the superconductor is forced to carry large 
currents. The possible role of impurities will be analyzed in a future paper.! 

Instead of using a locally self-consistent description, we have assumed that the gap is of 
the form (3) everywhere in the superconductor and zero in the normal lead. This assump- 
tion simplifies the numerics considerably and underlines the main effect of self- consistency, 
which is the introduction of a nonzero superfluid velocity v s = hq/m in the current- carrying 
configurations. In our scattering approach, we assume that quasiparticles are sent through 
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the incoming channels of the semi-infinite leads with chemical potentials that differ by eV. 
Then, we look for self-consistent values of |A| and q satisfying Eqs. (1-3) and asymptotic 
current conservation. We wish to emphasize that the occupations {f n } are those of a nonequi- 
librium population. The existence of a nonzero v s has the major implication that the energy 
thresholds for propagation across the superconductor are modified, as shown in the insets of 
Fig. lb. The minimum energy to excite a quasiparticle in the positive (+kp) and negative 
(— &f) branches becomes A + and A_, respectively, with A-t = \A\±h~VFq- |A| can in turn 
be supressed by current, in a manner that is very sensitive to the details of the scattering 
process. 

At zero temperature, electrons with energy < E < eV are injected from the normal 
side onto the NS interface. At low voltages (eV < A_), normal and Andreev reflection (AR) 
are the only scattering mechanisms. In a AR process, a hole is reflected on the N side 
and two electrons are transmitted to the superconductor forming a Cooper pair. Thus, one 
AR event contributes a charge 2e that is transported by the condensate. As the voltage is 
raised, q increases and, when eV > A_, Andreev transmission (AT) becomes also possible: 
one electron is transmitted to the superconductor, forming a Cooper pair and creating a 
quasi-hole. The amplitude for AT is proportional to r*t, where r (t) is the one-electron 
reflection (transmission) coefficient .0 One AT event contributes a net charge of e to the 
current, 2e going to the condensate and — e to the quasiparticle part. Thus, AT is less 
efficient than AR in that, for the same contribution to the total current, it displaces the 
condensate twice as much and thus costs more free energy. We remark that there is a range 
of voltages in which AT is the only quasiparticle transmission mechanism. By contrast, in 
a non self-consistent description, the AT and normal transmission (NT) channels have the 
same voltage threshold.0 

As V and q continue to increase, the threshold A_ becomes negative. In some cases, 
and especially at low temperatures, a type of gapless superconductivity (GS) with peculiar 
properties may survive. As shown in the right inset of Fig. lb, the quasiparticles that are 
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hidden below the e = line in the negative branch reemerge with positive energy in the 
opposite branch. For each quasiparticle n that satisfies Eq. (1), there is another solution n! 
with energy e n > = —e n and wave components u n i = —v* and v n > = w^.i Creating n', a is equiv- 
alent to destroying n, —a (a is the spin), and one is free to employ the n or n' description.^ 
The conventional choice is that which makes the energy positive, so that quasiparticles are 
rare at low temperatures. Following this convention, a quasiparticle n that acquires neg- 
ative energy when GS is reached must pass to be described by its n' counterpart. In the 
GS regime, low energy electrons coming from the N side can be normally transmitted as a 
quasi-electron in the lower positive branch. Being of the n'-type, this quasiparticle state has 
some unusual properties. Since u n v* n = — u n >v */, it is clear from Eq. (2) that the presence of 
n'-type quasiparticles tends to cancel the effect of conventional (n-type) quasiparticles. In a 
standard scenario without GS, the occupation of a quasiparticle state tends to decrease the 
value of the gap because the factor (1 — 2f n ) changes sign and tends to cancel the contribu- 
tion from unoccupied states. In the presence of GS, the group of n'-type quasiparticles acts 
in the opposite way. In front of already existing n-type quasiparticles, these unconventional 
quasiparticle states tend to depress the gap if they are empty and contribute to reinforce 
superconductivity if they are occupied. 

In the spirit of Ref.0, we have modeled the NS interface with a step function gap and 
a one-electron potential barrier Vq5(x) of effective strength Z = Vq/Tivf- The self-consistent 
current-voltage characteristics is computed for different values of Z. Except for a final 
prediction, we focus on the zero temperature case. We remark that the self-consistent 
solutions with |A| 7^0 have been checked to have lower free energy than the trivial solutions 
with zero q and |A|. We have taken a bandwidth of Ep = 5 eV, a cutoff energy of huJn = 0.1 
eV for Eq. (2), and a zero temperature, zero current gap of A = 1 meV. This yields a 
critical temperature of T c = 6.6 K. In Fig. la, we show the self-consistent I-V characteristics 
for an NS interface (e>0). Non self-consistent results0 are also shown for comparison. 

A detailed understanding of the physics behind the I-V curves shown in Fig. la would 
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require information on the dependence of | A|, v s , and the quasiparticle current on the applied 
voltage, which will be presented elsewhere.il Here, we only wish to note that, as the voltage 
increases, the first rise in the current is due to the opening of the AT channel (except in 
the Z = case, where it is forbidden because r = 0). The sharp decrease in / for Z = 
and 0.5 is due to the onset of the GS regime, where normal quasiparticle transmission into 
the unconventional branch begins to dominate, bringing down the value of the conductance 
from close to 2 (in units of2e 2 /h), typical of a transmissive NS interface, to near unity, as 
corresponds to reflectionless normal transport. This rather spectacular feature is however 
not likely to be measurable in practice, since it would be destroyed by the presence of 
a second normal lead. For example, in the transparent case, this second interface must 
send the value of the low V differential conductance back to its natural limit of 1. On 
the other hand, nonequilibrium effects that depend essentially on the total depopulation of 
incoming channels from the S side can hardly be manifest when a second normal lead sends 
quasiparticles from the opposite side. These considerations motivate our study of the NSN 
structure, for which we assume identical barriers and incoherent multiple scattering. By 
symmetry, the potential difference will be V/2 at each interface, and the flux of incoming 
electrons from the left N lead will be mirrored by the flux of incoming holes from the right 
N lead. 

In Fig. lb, we show the resulting I(V) curve for a NSN structure. Non self-consistent 
results obtained with analogous scattering assumptions are also shown. Comparison with 
Fig. la reveals indeed major changes. If Z = 0, the presence of a superconducting segment 
does not affect the conductance, since it adds nothing to perfect transmission. Being r = 0, 
AT is inhibited. When v = eV/A Q = 2, the system reaches GS, but, remarkably, |A| can 
be shown to remain unaffected.il This happens because, when the unconventional branch 
emerges, its two states are filled with very high probability and thus do not contribute to 
depress the gap. This surprising behavior contrasts markedly with that expected for situa- 
tions in which quasiparticles are in equilibrium. In such cases, |A| is predicted to decrease 
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very quickly as A_ becomes negativetm Here, a strong departure from equilibrium gives 
rise to an effective translational invariance that explains the small sensitivity of |A| to high 
voltages. 

For thicker barriers, we see that, in a non self-consistent calculation, the jumps in the 
current are much smoother and occur at higher V. For Z = 0.5, the AT channel opens for 
v ~ 1.3, causing a small decrease in |A|i, and GS sets in for v ~ 1.8. The onset of GS is 
almost inconsequential for the gap and total current behavior, for reasons similar to those 
of the transparent case. We conclude that, for NSN structures with low Z, the voltage 
is inefficient in destroying superconductivity, although this may remain quite marginal (in 
the GS regime, the actual current carried by the condensate is a small fraction of the total 
current). 

For Z = 1 and 2, the system bypasses the only- AT regime and jumps directly to GS for 
v ~1.7 and u~1.9, respectively. We attribute this behavior to the energetic cost of opening 
two AT channels. The superconductor finds more favourable to carry charge through the NT 
channel by going gapless. An important effect is that multiple normal reflection increases the 
residence time of the injected quasiparticles. At high Z, one can prove by invoking unitarity0 
that all the effective occupations within S tend to 1/2 (two of the four incident channels are 
empty), as opposed to the low Z case (see above) in which conventional (unconventional) 
channels are occupied with low (high) probability. This effect causes |A| to decrease with 
increasing Z, in marked contrast with the behavior found for the NS case.! 

The above discussion shows that, for the values of Z considered, superconductivity can 
exist at least up to voltages of order 5fcsT c / e. Internal thermalization of quasiparticles would 
tend to destroy this effect. In this sense, we can state that the survival of superconductivity 
at high voltages is a nonequilibrium effect. 

We would like to end this discussion with a few comments about the effect of temper- 
ature and dimensionality. Provided the superconductor can still be treated as quasi-one- 
dimensional (width much smaller than coherence and penetration lengths), the effect of many 
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channels is that of requiring higher q's to observe similar effects, since high-lying channels 
have a smaller effective longitudinal Fermi velocity. Also, the presence of many channels 
effectively breaks the translational invariance we found for very transmissive barriers.! As 
a consequence, the survival of superconductivity at high voltages is considerably weakened, 
in accordance with the behavior found in semiclassical nonequilibrium superconductivity.0 
Apart from smoothening the structure of the I(V) curves, the main effect of temperature 
is that of bringing down the scale of gap energies. With smaller gaps, the condensate needs 
higher values of q and the result is that self-consistency effects are even more important. 
A major consequence is that, as the voltage increases, the thresholds for the novel regimes 
discussed above are anticipated. Fig. 2 shows that the voltage threshold for the onset of 
AT, Vat, which determines the first peak in the differential conductance, is a decreasing 
function of temperature, even when measured in units of the zero current gap A (T). One 
may also note that Vat increases with Z, due to a global reduction of v s . 

In conclusion, we have performed a self-consistent calculation of the nonlinear transport 
properties of NS and NSN structures. A major effect of self-consistency is that of forcing 
the current- carrying superconductor to have a nonzero superfluid velocity. This causes the 
various quasiparticle channels to enter into action at different voltages. The study of the 
NS interface has served to identify the scattering processes separately and has allowed us to 
understand in detail the properties of the more realistic NSN system. The gap amplitude 
is depressed by the flow of current in a way that is very sensitive to the details of the 
scattering problem. An important nonequilibrium effect is that, at low temperatures, a type 
of gapless superconductivity may survive up to voltages considerably greater than /c^T c /e, 
although this effect may be weakened in the presence of many channels. The onset of 
Andreev trasmission is signaled by a peak in the differential conductance whose position 
relative to the zero current gap has been predicted to decrease with temperature and with 
the transmittivity of the barrier. 

Note added. After this work had been completed, we learned about the content of Ref.0, 
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where a related NSN model has been studied. In Ref.1121, one interface is always transmissive 
and multiple phase-coherent scattering is assumed. In spite of these differences, the model 
studied by Martin and LambertEl displays some common robust features like the splitting 
of voltage thresholds and the depression of the gap. 
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FIGURES 

FIG. 1. Solid lines show the self-consistent I-V characteristics of a NS (a) and a NSN (b) 
structure, for four values of Z. Dotted lines show the results from a non self-consistent calculation 
(see Ref. 9). The dashed line in (a) gives the NN result. Insets: schematic quasiparticle dispersion 
relation, e(k), in the (a) normal and (b) superconducting leads. Right inset in (b) shows gapless 
superconductivity. Filled (empty) circles indicate electron- (hole-) like propagation. 

FIG. 2. Position of the first peak in the differential conductance (caused by the onset of 
Andreev transmission), in units of the zero current gap, plotted as function of temperature for 
Z = 0.5, 1, and 2. 
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